Independence of adsorption transition point from bending en- 
ergy of a semiflexible polymer chain: Exact results 

Pramod Kumar Mishra 

Department of Physics, DSB Campus, Kumaun University, Naini 
Tal-263 002 (Uttarakhand) India. 

Email: pkmishrabhu@gmail.com 

Abstract : We analyze adsorption desorption phase transition behaviour 
of a linear semiflexible polymer chain on a geometrical constraint (s) under 
good solvent condition in two dimensions. The constraint is stair shaped 
impenetrable surface, for two dimensional space, the surface is a line and 
the polymer chain is confined in between two (A & B) such surface. Lattice 
model of fully directed self avoiding walk is used to calculate exact expression 
of the partition function when the chain has attractive interaction with any 
or both the constraints. We have found in each of the case that adsorption 
transition point of the polymer chain is independent of bending energy of 
the chain for different values of spacing between the constraints. Our result 
differs from predicted behaviour for adsorption of the semiflexible polymer 
chain on a flat surface where adsorption transition point depends on bending 
energy of the chain. 

PACS: 05.70.Fh, 64.60 Ak, 05.50.+q, 68.18.Jk, 36.20.-r 

1 Introduction 

Biopolymers (DNA & proteins) are soft object and therefore such molecule 
can be easily squeezed into the spaces that are much smaller than natu- 
ral size of the molecule. For example, actin filaments in eukaryotic cell or 
protein encapsulated in Ecoli [1] are the examples of confined molecules 
that serves as the basis for understanding molecular processes occurring in 
the living cells. The conformational properties of single biomolecule have 
attracted considerable attention in recent years due to developments in the 
single molecule based experiments [2, 3]. Under confined geometrical con- 
dition, excluded volume effect and effect of geometrical constraint compete 
with entropy of the molecule. Therefore, geometrical constraints can modify 
conformational properties and adsorption desorption transition behaviour of 
the molecules. 



The behaviour of a linear flexible polymer molecule under good solvent 
condition, confined to different geometries have been studied for past few- 
years [4, 5, 6] and semiflexible molecule based theoretical studies under con- 
fined geometry also find considerable attention in recent years [7, 8] and also 
see, references quoted therein. For example, Whittington and his coworkers 
[4, 5, 6] used directed self avoiding walk model to study behaviour of a flexi- 
ble polymer chain confined between two parallel walls on a square lattice and 
calculated force diagram for surface interacting polymer chain. Rensburg et. 
al [6] showed through numerical studies using isotropic self avoiding walk 
model that force diagram obtained for surface interacting polymer chain 
confined in between two parallel plates have qualitatively similar phase dia- 
gram to that has been obtained by Brak et. al [4] for directed self avoiding 
walk model. 

However, in present investigation, we have considered a linear semiflex- 
ible polymer chain confined in between two one dimensional stair shaped 
impenetrable surface (geometrical constraint) under good solvent condition 
and we have discussed adsorption phase transition behaviour of an infinitely 
long linear semiflexible polymer chain on the constraint (s). Such investiga- 
tion may be useful to understand biological process of macromolecule and 
membrane interactions as well as such process of biotechnology as DNA in 
micro-arrays and electrophoresis. 

To analyze adsorption phase transition behaviour of the semiflexible 
chain, confined in between the constraints, we have chosen fully directed 
self avoiding walk model introduced by Privmann and his coworkers [9] and 
used generating function technique to solve the model analytically for dif- 
ferent values of spacing between the constraints. The results so obtained is 
used to discuss adsorption phase transition behaviour of the polymer chain 
on the stair shaped surface and also to compare the results obtained for 
adsorption transition of the polymer chain on a flat surface [10, 11]. 

Since, the constraint is an attractive surface, therefore,it contributes 
an energy e s (< 0) for each step of the fully directed self avoiding walk 
making on the constraint. This leads to an increased probability defined by 
a Boltzmann weight uj = exp(—e s /kBT) of making a step on the constraint 
(e s < or oj > 1, T is temperature and ks is the Boltzmann constant). The 
polymer chain gets adsorbed on the constraint at an appropriate value of oj 
or e s . Therefore, transition between adsorbed to desorbed phase is marked 
by a critical value of adsorption energy e s or oj c . The crossover exponent (4>) 
at the transition point is defined as, N s ~ N^, where N is the total number 
of monomers in the chain while N s is the number of monomers adsorbed on 
the surface. 

In this paper, we have analytically solved fully directed self avoiding 
walk model to investigate dependency of adsorption transition point of an 
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infinitely long linear semiflexible polymer chain when adsorption of the chain 
occurs on the constraint (s) and used the results so obtained to compare with 
the case when adsorption of a semiflexible polymer chain occurs on a flat 
surface [10, 11]. 

The paper is organized as follows: In Sec. 2, lattice model of fully 
directed self avoiding walk is described for a linear semiflexible polymer 
chain confined in between the constraints for a particular value of spacing 
between the constraints under good solvent condition on a square lattice. 
In sub-section 2.1, we have discussed adsorption behaviour of the polymer 
chain when constraint A is having attractive interaction with the semiflex- 
ible polymer chain. Sub-section 2.2 is devoted to discuss the adsorption of 
the semiflexible polymer chain on the constraint B. While, in sub-section 
2.3, expression of the partition function of the polymer chain is obtained 
for the case when the chain is having attractive interaction with the both 
the constraints. Finally, in Sec. 3 we summarize and discuss the results 
obtained. 

2 Model and method 

Lattice model of fully directed self-avoiding walk [9] on a square lattice 
has been used to analyze dependency of adsorption transition point of an 
infinitely long linear semiflexible homopolymer chain on its bending energy 
when the chain is confined in between two impenetrable stair shaped surface 
under good solvent condition (as shown schematically in figure 1). The 
directed walk model is restrictive in the sense that the angle of bending has 
unique value, that is 90° for a square lattice and directedness of the walk 
amounts to certain degree of stiffness in the walks of the chain because all 
directions of the space are not treated equally. Since, directed self avoiding 
walk model can be solved analytically and therefore it gives exact value 
of the adsorption transition point of the polymer chain. We consider fully 
directed self avoiding walk (FDSAW) model, therefore, walker is allowed 
to take steps along +x, and +y directions on a square lattice in between the 
constraints. The walks of the chain starts from a point O, located on an 
impenetrable surface and walker moves through out the space in between the 
two surface (as we have shown schematically in figure (1) that a walk of the 
polymer chain confined in between two surface for one value of separation 
n(= 3) between the two surface or the constraints). 

The stiffness of the chain is accounted by associating a Boltzmann weight 
with bending energy for each turn in the walk of the polymer chain. The 
stiffness weight is fc(= exp(—(3eb); where j3 = is inverse of the temper- 
ature, eb(> 0) is the energy associated with each bend in the walk of the 
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Figure 1: This figure shows a walk of an infinitely long linear semiflexi- 
ble polymer chain confined in between two constraints (impenetrable stair 
shaped surface). All walks of the chain starts from a point O on the con- 
straint. We have shown three different cases viz. (i), (ii) and (iii) having 
separation (n) between the constraints along an axis three monomers (steps). 
The separation between the constraints have been defined on the basis of 
the fact that how many maximum number of steps a walker can successively 
move along any of the +x or +y direction. In the case l(i), the constraint A 
is having attractive interaction with the monomers of the chain, in 1 (ii) only 
constraint B is having attractive interaction with the monomers of the chain 
while in 1 (iii) both constraints are shown to have attractive interaction with 
the monomer of the polymer chain. 
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chain, kb is Boltzmann constant and T is temperature) . For k = 1 or e& = 
the chain is said to be flexible and for < k < 1 or < et < oo the polymer 
chain is said to be semiflexible. However, when q, — > oo or k — > 0, the chain 
has shape like a rigid rod. 

The partition function of a surface interacting semiflexible polymer chain 
can be written as, 

a// walks of N steps 

where, iVj, is the total number of bends in a walk of N steps (monomers), 
iV s is number of monomers in a N step walk, lying on the surface, g is the 
step fugacity of each monomer of the chain and oj is Boltzmann weight of 
monomer-surface attraction energy. 

2.1 Adsorption of the semiflexible polymer chain on the con- 
straint A 

The partition function of an infinitely long linear semiflexible polymer chain 
confined in between the constraints (as shown schematically in figure 1-i) 
and having attractive interaction with the constraint A can be calculated 
using the method of generating function technique. The components (as 
shown in figure 2) of the partition function, Z A (k,0Ji) (we have used here 
suffix three because in figure 1-i case maximum step that a walker can move 
successively in one particular direction is three and lo± is Boltzmann weight 
of attraction energy between monomers and the constraint A) of the chain 
can be written as, 

X A = si + k Sl Y 3 A (2) 

where, s\ = uj±g. 

Xf = g + g(X? + kY 2 A ) (3) 
x£ = g + g(X A + kY 1 A ) (4) 
Y A = g + kgX A (5) 
Y A = g + g(kX A +Y A ) (6) 

and 

Y A = s 1 + s 1 (kX A + Y A ) (7) 

On solving Eqs. (2-7), we get the expression for X A (k, u>i) and Y A (k, oo\). 
In obtaining the expression for X A {k,uj\) and Y A (k,uj\), we have solved a 
matrix of 2nX2n (n = 3, for present case i. e. figure 1-i). Thus, we 
have exact expression of the partition function for an infinitely long linear 
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Figure 2: The components of the partition function is shown graphically 
in this figure. Term X A (3 < m < n) indicates sum of Boltzmann weight 
of all the walks having first step along +x direction and suffix n indicates 
maximum number of steps that walker can successively take along +x direc- 
tion. Similarly, we have defined where first step of the walker is along 
+y direction. In this figure, (ii) and (iii) graphically represents recursion 
relation for Eqs.(2&6) respectively. 




semiflexible polymer chain confined between the constraints and having at- 
tractive interaction with the constraint A (as shown in figure I-i) is written 

as, 

rrAii \ ~v A f i \ t^4/i \ u\ + u 2 + u 3 + 2k A s\x A - 2k b s\x A 
Z^(k, Ul ) = Xf(k, Ul ) + Y 2 A (k, Ul ) = _ l + k , s 2 x l +UA — 

(8) 

where, 

u\ = s± — ks\ — x — ks\x + k 2 s\x — x 2 — kx 2 — ks\x 2 + 2k 2 s\x 2 — ks 2 x 2 
u 2 = -k 2 s\x 2 + 3k 3 sjx 2 - kx 3 + k 2 x 3 - k 2 sjx 3 + 2k 3 six 3 - k A s 2 x 3 
u 3 = —kx A + k 3 x A — ks\x A + k 2 six A + k 3 six A — k A six A — 2k 2 s 2 x A + 2k 3 s 2 x A 
and 

u 4 = -k A {x A + 2s 2 {x 2 + x A )) + k 2 (x 2 (2 + x 2 ) + sj(l + x 2 + x A )) 
From singularity of the partition function, Z^(k,ui), we obtain critical 
value of monomer-constraint A attraction energy, 

Ucl ~ yJk 2 g 2 +k 2 g' i ~2k' i g' i +k 2 gS-2k*gZ+kZgZ ' 

required for adsorption of an infinitely long linear semiflexible polymer 
chain on the constraint A. We have obtained the value of oj c \ = 1, when 
we substitute value of g c in the expression of uj c \ corresponding to all pos- 
sible values of k[= exp(— /3e&)] or bending energy e& for which an infinitely 
long linear semiflexible polymer chain can be polymerized in between the 
constraints. 
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2.2 Adsorption of the semiflexible polymer chain on the con- 
straint B 



The partition function of an infinitely long linear semiflexible polymer chain 
confined in between the constraints (as shown schematically in figure 1- 
ii) and having attractive interaction with the constraint B is calculated 
following the method discussed in above subsection. The components of the 
partition function, Z B {k,ui2) (where, L02 is Boltzmann weight of attractive 
interaction energy between the monomers of the chain and the constraint 
B) of the chain can be written as, 

X B =g + kgY B (9) 

X B =g + g(X B + kY 2 B ) (10) 
X B = s 2 + s 2 (X B + kY B ) (11) 



where, S2 = ^29- 



and 



Yf = s 2 + ks 2 X B (12) 
Y 2 B = g + g(kX B + Yf) (13) 



Y B = g + g(kX B + Y B ) (14) 

On solving Eqs. (9-14), we get the expression for X B (k, co 2 ) and Y 2 B (k, u^)- 
In obtaining the expression for X B (k,uj2) and Y B (k,oj 2 ), we have to solve 
a matrix of 2nX2n (n = 3, for figure 1-ii case). Thus, we obtained exact 
expression of the partition function for an infinitely long linear semiflexi- 
ble polymer chain confined between the constraints and having attractive 
interaction with the constraint B (as shown in figure 1-ii) is as, 

7%,n-y%nUV% f n- -x(g 2 (l + kx 2 - k 2 x 2 ) + 2u 5 + ks 2 2 u 6 ) 

Z 3 {k,U 2 ) - X 1 {k,LJ 2 ) + Y 2 {k,UJ 2 ) _ 1 + k 6 s 2 x 4 + Uj 

2 ' (15) 

where, 

u 5 = 1 + k 2 (-l + x)x 2 - k 3 x 3 + kx{\ + x) 

u 6 = i + x + x 2-2k 3 {-l + x)x 2 + 2k 4 x 3 -k 2 x{2 + 3x + 3x 2 ) + 2k{-l + x 3 ) 
u 7 = -k 4 (x A + 2s 2 2 (x 2 + x 4 )) + k 2 {x 2 {2 + x 2 ) + s 2 {l + x 2 + x 4 )) 
From singularity of the partition function, Z B (k,oj 2 ), we could obtain 
critical value of monomer-constraint B attraction energy, 

Uc2 ~ ^V+^^^V+fcV^V+P^e _ Wcl ' 

required for adsorption of an infinitely long linear semiflexible polymer 
chain on the constraint B. In this case too, we have value of lo C 2 = 1, for 
all possible values of bending energy or stiffness of the semiflexible polymer 
chain. 
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2.3 Adsorption of the semiflexible polymer chain on the con- 
straints ALB 



The partition function of an infinitely long linear semiflexible polymer chain 
confined in between the constraints (as shown schematically in figure 1- 
iii) and having attractive interaction with both the constraints (A&iB) is 
calculated following the method discussed above. The components of the 
partition function, Z% (k, 003,1^4) of the chain can be written as, 



where, S3 = u 3 g. 



here, S4 = u^g. 



and 



Xf = s 3 + ks 3 Y 3 c (16) 

xC = g + g(xC + kY 2 c ) (17) 

Xf = s 4 + s A (X% + kYf) (18) 

Yf = Si + fcs 4 Xp (19) 

Y 2 C = g + g{kX% + Yf) (20) 



Yf = s 3 + s 3 (kX^ + Yf) (21) 

On solving Eqs. (16-21), we get the expression for X^ (k, uj 3 , uj 3 ) and 
Yf(k, uj 3 , U4). In obtaining the expression for Xi(k, ^3,(^4) and Y 2 c (k, ^3,(^4), 
we have solved a matrix of 2nX2n (n = 3, for figure 1-iii case). Thus, we 
have exact expression of the partition function for an infinitely long linear 
semiflexible polymer chain confined between the constraints and having at- 
tractive interaction with the constraints (as shown in figure 1-iii) is written 
as, 

Z 3 (k,^) = X, (k,u 3 ,u 4 )+Y 2 (fc, 0^,0,4) = - _ 1 + k e s 2 s 2 x 2 + Ul2 + Ul3 

(22) 

where, 

us = 1 + S4 + kx — (—1 + k)ks\{\ + x + kx) 

Ug = 1 - k 2 s\x 2 + k A s\x 2 + k(l + s\)x 2 - k 2 {x 2 s\{\ + x 2 )) 

nio = ks 2 (l + x + S4xk 2 s 2 (l — 2x)x + 2k 4 s 2 x 2 + kx(— 1 — S4 + 2x)) 

Ull = ks 2 3 {kxs\{\ + 2x) - k 2 {2x 2 + s\{\ + 2x + 2x 2 ))) 

Ul2 = k 2 (x 2 + s\{\ + x 2 ) + s 2 (l + (1 + sl)x 2 )) 

and 

«13 = -k\s 2 x 2 + s 2 {x 2 + s\(l + 2x 2 ))) 
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Prom singularity of the partition function, Z% (k, L03, 0J4), we obtain crit- 
ical value of monomer-constraints attraction energy, 

Jl - k 2 g 2 - k 2 g 2 u 2 - k 2 g 4 coj + k 4 g 4 uj 
w c3 = -, =, (23) 

^k 2 g 2 + k 2 g 4 - k 4 g 4 - k 4 g^uj\ - 2k 4 g & oj 2 + k 2 g^oj\ + k^g^ojj 

required for adsorption of an infinitely long linear semiflexible polymer 
chain on the constraints A when both the constraints have attractive inter- 
action with the chain. 

On substitution of following value of uja in Eq. (23) to get value of 
W C 3 = 1. 

The method discussed above can be used for different values of n and 
size of the matrix needed to solve in calculating partition function of the 
chain confined in between the constraints is 2nX2n. We have calculated 
exact expression of the partition function for n (3 < n < 19). 

We have found that adsorption transition point of an infinitely long 
linear semiflexible polymer chain on the constraint A, B and simultaneously 
on both the constraints A&iB has value unity. This fact is true for all the 
chosen values of k or bending energy (and 3 < n < 00) for which an infinitely 
long polymer chain can be polymerized in between the constraints. 



3 Result and discussion 

We have considered an infinitely long linear semiflexible homopolymer chain 
confined in between two impenetrable stair shaped surface (constraint) on 
square lattice under good solvent condition. We have used fully directed 
self avoiding walk model to study adsorption phase transition behaviour of 
the polymer chain on any of the two constraints {A&i B) and simultaneous 
adsorption of the polymer chain on both the constraints AhB. The gener- 
ating function technique is used to solve the model analytically and exact 
expression of the partition function of surface interacting semiflexible poly- 
mer chain is obtained for different values of spacing (3 < n < 19) between 
the constraints. 

It has been found that adsorption of the polymer chain occurs on the 
constraints at a value uj c \ = uj C 2 = wc3 = 1 for all possible values of k or 
bending energy of the chain for which an infinitely long linear semiflexible 
polymer chain can be polymerized in between the constraints. The critical 
value of uj is unity for all the considered cases by us for different values 
of spacing between the constraints (3 < n < 19). This result is obvious 
because walks of the chain is directed along the constraint (s), therefore, 
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partition function of the chain is dominated by walks lying on the constraints 
and adsorption of the chain will occur even there is no attraction between 
constraint (s) and monomer of the chain. 

However, in the case of adsorption of an infinitely long linear semiflexible 
polymer chain on a flat surface, the adsorption transition point is found to 
depend on the bending energy or stiffness of the chain. In this case, stiffer 
chains adsorption occurs at a smaller value of monomer surface attraction 
than the flexible polymer chain [10, 11]. 
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